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1. Introduction 



Experimental knowledge of the flow field generated by 
rotating turboimpellers is essential for the research and 
development of turbomachinery. This information is used to 
refine design methods, develop new flow models which include 
secondary flow and tip clearance effects, and especially to 
verify computer programs designed to calculate flow throuoh 
rotating blade rows. 

Laser velocimeters have been used successfully in recent 
years to measure the flow inside and downstream of rotors 
(see Ref. 1) . Certain disadvantages have become apparent, 
however. The laser techniques are reliable only in the hands 
of experienced investigators, the pressure field remains 
unknown, and usually the measurement of more than two compo- 
nents of the velocity field is complicated and expensive. 
Furthermore, it is difficult to perform measurements close to 
walls. Development of alternative techniques to overcome 
these deficiencies, as well as to achieve redundancy in mea- 
suring the flow field, are reasonable and worthwhile tasks. 

This report describes a particular method and the com- 
putational support necessary to measure the flow field behind 
an impeller in the stationary, bladeless gap. 
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2 . 



Description of Method 



The following method requires two semiconductor pressure 
probes along with a technique for synchronized sampling for 
determining the fluid velocity vector downstream of a rotor. 

The two probes (see Fig. 1) are positioned inside the 
machine casing so they will, in turn, intercept periodically 
the same part of the flow leaving a particular passing rotor 
passage. Each probe reading is sampled when the designated 
blade passage reaches a desired position relative to the 
probe. Synchronization is achieved through a suitable method 
(Ref. 2,3). 

Four quantities are needed to determine the velocity 
vector: yaw angle, pitch angle, static pressure and total 

pressure. Accordingly, four measurements must be made to 
evaluate these unknowns. By rotating the probes about their 
tips, pressure readings in four different directions can be 
taken, and the data used to calculate the velocity vector. 
Computer program VELOCITY, given in Appendix II, was developed 
to perform the somewhat arduous calculations. 

The geometries of the two probes are shown in Fig. 1. 
Before being used, the probes must be calibrated so their 
responses to flows coming from different directions are known. 
A highly directional probe is desired to increase the accuracy 
in finding the yaw and pitch angles, and consequently the 
velocity magnitude. The following method is recommended for 
calibrating each probe - 
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PROBE TYPE A 



PROBE TYPE B 



Figure 1. A type and B type probes 



3 



1 . 



Establish a steady, controlled flow of fluid, and 
determine the velocity vector at a certain region of 
the flow. 

2. Position a probe in the flow and rotate the tip so 
that a sequence of pressure readings are taken for 
a constant yaw angle and a varying pitch angle. 
Repeat the procedure at a new yaw angle using the 
same pitch angles. The result will be an array of 
pressure readings corresponding to a set grid of 
yaw and pitch angles (Fig. 2) . 

3. From the known flow velocity and pressure readings, 
a coefficient of pressure can be calculated for each 
angle set: 



p-p 

C = — 

P Prr-P. 



where : 



Cp = Coefficient of pressure 
p = pressure reading 
P s = static pressure of flow 
p^p = total pressure of flow 



The table of C^'s as we H as the yaw and pitch angles which 
correspond to them are now in the form required for input to 
program VELOCITY. 

The probe calibrations should be insensitive to Mach 
number and pressure, and are not valid for supersonic flows. 
Should any significant variations in C^ be observed for 
different flow conditions, further calibrations will be re- 
quired and an additional iteration scheme added to the com- 
puter program. 
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Experience with the two-probe technique has shown that 
excellent results are achieved when a probe type A is rotated 
to the three positions +25°, 0°, -25° yaw.at 0° pitch) ,and 
probe type B is used at 0° yaw and 25° pitch, Fig. 3) . 

The two-probe technique is strictly applicable only to 
periodic flows. However, data obtained on successive rotations 
of the rotor can be averaged to eliminate non-periodic fluctua- 
tions. This was effective for tests reported in Ref. 2., where 
a single probe was used to establish the peripheral blade- to- 
blade distribution of flow yaw angle. 

It is noted that the method reported here is a further 
development of that reported earlier in Ref. 6, and overcomes 
some of the earlier limitations. 
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3. Theory 



The velocity vector for a three-dimensional flow can be 
described with three scalar quantities. The nature of the 
problem suggests using two angles (a yaw angle and a pitch 
angle) , and the magnitude of the velocity (Fig. 4) . 

Since pressures and not the velocity are measured, the 
static and total pressures must first be determined, and 
Eq. (1) used to evaluate the velocity. 



P T 

p s 



(1 + 



111 m 2jY/Y-1 



( 1 ) 



Altogether, four unknowns need to be evaluated: the yaw and 

pitch angles, and the total and static pressures. 

Four equations are needed to determine the four unknowns. 
They are derived from the four pressure readings, each pres- 
sure reading having been taken in a different direction as 
described above. The following equations for the coefficient 
of pressure can be written: 




p r p s 

P T _P s 



i = 1. . 4 



( 2 ) 



The C .'s are a function of the orientation of the probe 

pi * 

relative to the flow; i.e. , for a given flow the measured C 's 

P 

will vary measureably as the probe is turned into and away 
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Figure 4. Velocity Vector V 
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from the flow. Each "probe"* will have its own C charac- 

ir 

teristics determined experimentally. The result will be a 
table of Cp vs. yaw and pitch angles for each probe. 



C pi “ fUnCti ° n <“ Ri ' W 



i = 1. .4 



(3) 



For realistic problems, only one point (a , <j>) exists 

where the C . 's in Eq. (2) will equal the C . , of Eq. (3) for 
pi ^ ^ pi s 

the four probes' pressure readings. 

The probes' characteristics (C^'s) are in tabular form 
because they cannot be represented analytically due to the sten 
effect and production inaccuracies. Therefore, a numerical 
solution to the problem is required. The procedure chosen for 
solving the problem is a systematic trial- and-error search 
process, essentially a convergence scheme on two variables: 
yaw angle and pitch angle. 

The flow direction is assumed to fall within some set of 
bounds, defining the search area for yaw and pitch (Fig. 5). 

By setting up a grid of points in this region and checking how 
well each point satisfies the criteria of equality of coeffi- 
cients of pressure (C .’s) calculated with Eqs. (2) and (3), 
the point with the smallest error can be found and used as a 
first approximation to the solution. Repeating this procedure, 
only with a smaller grid and search region, will result in a 
better approximation. This sequence, represented in Figs. 6 



* Here, the term "probe" refers to a particular probe type in 
a particular position. 
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and 7 is repeated until either the desired accuracy is reached 
or fatigue sets in. Program VELOCITY, described in the follow- 
ing section, was written to perform these calculations. 
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Figure 6. Illustration of the Search Procedure 
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Figure 7. Flow Chart of the Search Procedure 
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4. Program VELOCITY 



Program VELOCITY was written to perform the calculations 
outlined in the previous section. A description of the program 
and its subroutines is given below. Fig. 8 summarizes the 
major sections and organization of the program. 

For each run, program VELOCITY reads the calibration 
tables for the two probes from files outside the program. 

(Input formatting is discussed in Appendix V.) Subroutine 
INPUT performs the necessary work, and can be modified to 
accommodate different input schemes if desired. 

The fluid temperature and molecular weight are entered 
next. These properties are assumed to remain constant through- 
out the run. 

The settings for each pressure reading are read next. A 
setting contains the following data: probe type (A or B) , yaw 

angle setting, and pitch angle setting. Again, these settings 
will not change for the duration of the run. 

Finally, the four pressure readings are entered. 

The first scan is initiated and covers the entire region 
of expected flow directions, -40° to + 40° in both yaw and 
pitch angles in the present case. Points are chosen every 5°, 
each point representing a unique pair of yaw and pitch angles. 
For each point, a static pressure, a dynamic pressure, and an 
error are calculated by the scheme described below. 

A point, say (a , 4> ) is tested; i.e., a test is performed 
to prove whether assumed flow, oriented a degrees yaw and 
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Figure 8. 



Flow Chart of Program VELOCITY 
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degrees pitch relative to the laboratory reference frame, 
corresponds to the four pressure readings. The direction of 
the flow relative to each probe setting is calculated. For 
probe setting i , oriented at (a^ , <f>^) relative to the 
laboratory, the assumed flow approaches at a relative angle 
of : 



4 > 



Ri 






(5) 



where (ct D . , <f> .) are the yaw and pitch angles respectively 
of the assumed flow relative to probe setting i . The 
calibration table for the probe used in setting i is con- 
sulted and a C (a . , returned. Subroutine CPCAL locates 

p Kl Kl 

or calculates the desired C values in the table. The scheme 

P 

used in CPCAL is a search technique to find the values of yaw 
and pitch surrounding the desired point, and then a linear 
interpolation over these four points as shown in Fig. 9. 

Eq. (2) can be rewritten in the form 

(C pi >P T + (1 -V P s = P i i=1 -' 4 (6) 

the only unknowns being p T and p . With four equations and 
two unknowns, the problem will be inconsistent unless the true 
a and <f> were chosen. Accordingly, the following schemes 
were used to evaluate p g , p T and an error. 
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Figure 9. L i i . ea ^ interpolation between four points to find 



Define : 



C 

_R 



4 

E C n 

i=l p i 






C = minimum of (C , C , C , C ) 
P m P 1 p 2 p 3 p 4 



corresponding to the C chosen above. 



m 



(c p ,p T + <4-C p) p s = P 



and also 



(C )p^ + (1-C )p = p 
P™ T P m s ^m 

m m 



These two equations can be solved for p T and p £ 



p ( l-c ) - 



P T = 



^m 

“C - 4c 
P P 



p ( 4-C ) 

p 



m 



C (Pj - C (p) 



m 



P s = 



p m 



C - 4C 
P P 



m 



4 p . -p 

Error = L ABS(C - - - ~ ) /4 

i=l P i p T~ P s 



(7) 



( 8 ) 



(9) 



( 10 ) 



(ID 



( 12 ) 



(13) 



19 



These schemes were chosen for two reasons: 



1) They used all the available data to derive an error 
which would effectively represent the accuracy of 
the guess. 

2) No singularities in the calculations can occur except 

for the case of four equal C 's (which physically 

P 

represents trying to find an intersection point 
among four parallel lines) . If the measurements are 
taken in the suggested directions, this anomalous 
point will not appear. 

For each point guessed in the initial scan, an error is 
calculated and the point with the smallest error is saved. A 
new, finer search grid is composed using this point as the new 
origin. The boundaries of the new grid are the points from 
the old grid which were closest to this new origin. Referring 
to Figure 10, if x represents the true solution, the new 
boundary would be formed by the points marked B-I, and the new 
grid-width would be one third as large. This factor was chosen 
to minimize the number of guess evaluations. (The first scan 
contains a large number of guesses in order to correctly isolate 
the general region of the solution) . 

The search procedure is performed on each new grid, and 
the process repeated until the grid width is less that 0.5°. 
After the final scan, the best guess is used to calculate the 
flow velocity and Mach number. The results are printed out 
and the next four pressures requested. If no values are 
entered (end of data set) , the program ends. 
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• x TRUE SOLUTION 

Figure 10 . Defining new grid boundaries from the nearest neighbors 
of the point with the smallest error 
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5. 



Discussion 



Extensive tests with program VELOCITY have led to the 
observations and suggestions listed below: 

1. Excellent results are achieved when the probe settings 

are at (yaw, pitch) angles of (-25,0) , (0,0) , (25,0) 

and (0,25) degrees. This corresponds to a rotation 

of probe type A from -25° to 0° to 25°, and one read- 
ing from probe type B at (0,25). Poor results were 
achieved for the symmetric case of readings at (+25,0) 
and 0,+25) degrees. 

2. Highly directional probes increase the accuracy of 
the procedure, especially if the variation is 
significant when the flow is nearly head-on To 
achieve these characteristics, the following design 
suggestion is offered. The probe can be formed with 
a spherical tip, the pressure tap being located in 
the center. To prevent damage to the sensitive trans- 
ducer located behind the pressure tap and to improve 
the frequency response, the void between the pressure 
tap face and the transducer should be filled with an 
an appropriate liquid and the opening of the pressure 
tap sealed with a thin, low-inertia membrane. 

3. Higher accuracy naturally results if more calibration 
points are taken for the probes' tables. The linear 
interpolation scheme can be replaced by the second 
order scheme offered in Appendix 5 (if no significant 
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anomalies occur in the calibrations) , the second order 
method requiring fewer calibration points (say every 
15°) than the linear method (every 5° or 10°) . 

4. The use of two probes of relatively simple geometry 
in periodic flow is less cumbersome and complex than 
the use of five-hole probes (Ref. 4). 



23 



Notation Summary 



c 

p 




Coefficient of pressure 

C is a function of a 
P 


and 


c)> , C = C (a , 
P P 


C pi 


- 


Coefficient of pressure 
C pi - C p^ a Ri ' ^Ri^ 


for 


probe setting i 


fp 


- 


Sum of 


the four C . 1 s 
Pi 






c p m 


- 


Minimum 


of the four C . 

P 1 


’s 




c pn 


- 


c p ( a i ' 


V 






Cp 12 


- 


c P (a i ’ 


V 






c 

p 21 


- 


C p ,a 2 ' 


<£-]_) 






c 

P 22 


- 


C p ( “2 ' 


4> 2 ) 






P 


- 


Pressure (all pressures 


are 


absolute) 


P. 

1 


- 


Pressure read from probe setting i 


P s 


- 


Static 


pressure 







Total pressure (stagnation pressure) 



P 

P 



m 



Sum of the four pressures (Pi's) 



Pressure at the setting where C occurred U.e 

P 

^m 

P = P. , where i = m , defined in C ) 
mi p 



V 


Velocity magnitude of the fluid particle 


V 


Fluid velocity vector 


a , <f> - 


Yaw, Pitch angles 


a i ' *i " 


Yaw, Pitch angles for probe setting i 


a Ri ' ^ Ri 


- Yaw, Pitch angles for the assumed flow direction 
direction relative to the probe setting 


P 


fluid density 
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ESTAELISH SCANMNC RANGE, GRID WIDTH, AND INITIALISE CRRCR 
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COO 000 I 0 
C CO 00020 
C CO 01 0.3 0 
CC0G1040 

: c o o o o 5 o 
cocjoe: 
ocoojon 
' oo cocao 
: c o c o c i c 

n OOCO 100 

•: o i c o i o 
: c o c o i 2 c 
: C) co i 30 
c : i co i a o 
:coc o iso 
c o o : o 1 1 o 
r c i o o i 7 i 

U C 3 C C 1 * 0 
: : j c 1 1 oo 
: oi co2^c 
; o o c o ? i o 
:cooo 2 ?o 

CC1C0230 
31000240 
C OO C0250 

: c ■) o o 2 ft j 

0 G 0 C 0 2 7 C 
00701280 

:cooo 2 oo 

70100 ? 30 
C C 0 C C 2 1 C 
C OQC 0 320 
0 00 CO 3 30 

: oo C034C 

C01C0350 
COO CO 360 
C0100370 
0 CO CO 380 
00000390 
OOCC0400 

: oo co4 i o 
: o o c o 4 2 o 

CCO C043C 
COO 00440 
CCOC1450 
C0000460 
C0000470 

cooco4ao 

CC0CC490 

coocosoo 

C000051 0 
CCOC0520 
00000530 
C1000540 
C00C0S5C 
00000560 
C00C0570 
CCOOO 57 1 
CIO 00572 
C0OC0573 
C 0 0 CO 574 
C 0000 575 
C 0000576 
00000530 
CIO C0590 
COO 00595 
C 00 0 0 5 9 6 
C CO C0597 
C OC C0596 
00000599 
COO C 0 6 0 0 

coo ooe i o 

C CO 00620 
0 10C064C 
00000650 
C 00 00660 
C0000670 
00000680 



nc r.c rr.rr fir.nncn < r> r,r nnnnr.f, nr>rnr>r. nnr.r.n 



VELOCITY 



AYI\, AM AX = V IN IMV , MAXIMUM 
I is • FMAX - viMf.LV, VAXIMJM 
'Jhl = CMC UCTh 

= MMM;^ ERRCR ECLNC 

4VIh=-/»C. 

A V A > = 4 C « 

FVlf^ = -4Q. 
p v A > = 4 o • 

DEl = 5 . 

150 tKTMN-lCCCCO. 

STA^T SCAN FTCCEoUHC 

x = ancle cues s 

v = PI TCP ANCLC C-LF3S 



Y A a ANGLES 
PITCH ANCLCS 



SC FAR 



Y - P N IN 
X - A Y I N 



EJV 


- STlEES 


TEE 


SUM CP THF 


FCLR CP* 5 READ 


FROM '3 Y CPCAL 


: o 3 coe )o 


« SUM 


- S T C R E S 


TEE 


SUM CF THE 


FCUR INPUT FRFS 


SURFS 


30 JC 0900 


N I N 


= S T C R E S 


TEC 


M IN I MUN CP 


VALUE FDR TEIS 


G L E S S 


COO CO 9 1 0 


N I N 


= STCEFS 


TEC 


PRESSURE CORRESPOND I NG TC 


TEC MIN IMUM CP 


J 03 009 1 5 














C 30 C 0520 


CFZ 


L v = c . 










: o o c o 5 i o 


PR 3 


S U v _ C . 










C 30 C 0940 


CPM 


IN -5 • 










0 C 3 C 0 9 5 0 



I to 



1 70 



ST ART TK ANALYSIS BY TINOING THE CP VALUES FRCM TEC T4ELC CPCAL# 

A.MC C/ALLATING CPSLM, CP«IN. AND PRSSUM 

DC 2 0 r K = 1 , 4 

> rT ->-ALP(K ) 

Y 9 = Y- PH I ( K ) 

I K ( NPr F ( K ) • E C • 1 ) CALL CPCAL ( N ALPH 1 .NPh I 1 • PR C i . CP 1 .XP.YP.CP(K) • I F L ) 

I E ( NPRC < H ) . TO . 2 ) CALL CPC AL<N A L PH 2#NPHI2#PRP2,CF2.XR,YR,CP(K),IFL) 
IE(IFL.NF*C) CCTC 2 SO 
CF5LV=CF C LVKP(K) 

PHDSUV=FRSSC^+FRESS(K) 

I F ( Cl- V I N . t T .CP ( K ) ) GOTO 200 
CPV IN-CP ( K ) 

REN 1N-PFES S (K ) 

ZOO CCMINUE 

FkCV TEE AECVZ DATA, CALCULATE A TOTAL AND STATIC FFESSURC 

PTI = A CHAR AC TCfxi ST IC TCTAL PRESSURE E C R TEIS YA*, PITCH 

Pb S - A lH A F ALTER 1ST IC STATIC PRESSURE E CR THIS YA*, PITCH 



COOCOCRC 
00000700 
C 3 0 0 0 7 1 0 
000 03 720 
JOOOO 730 
0 00 CO 740 
COO 3 0 753 

: o ) c o 7 6 o 

030C977C 
.’0 3 00 7H3 
)0 3 00 790 
0 3CC0900 
C 3 0 C C P 1 C 

c o o c o e 2 o 
c o j c c e 3 o 

C 0030940 
OOOCOdSO 
: 00 C 095 G 
C 0 3 0 J 5 T 0 
■ ’ 0 0 C 0 5 9 0 



3F_N CM =CFS L N-4 . *CP VI N 

PTT = ( PRS SUN* < 1 #-CPM I N ) - PPMIN*< 3.-CRSLW) 

•>LS-( CPSt^PFMN - PES SUM*CPM IN) / OENCV 



) /CENCM 



CALCULATE a CHARACTERISTIC ERRCR and COMPARE fc I TF T E* C 
PREVICLSI.Y P C L NO SMALLEST ERRCR 

1F( FTT.LE.F5S) CCTC 250 
b»u R = C • 

DC r 25 I E - 1 « 4 

!25 ZEE F = L R‘ R E f AfcS(CP(IE) - ( PFESS ( I R ) ~ PS S ) / { FT T-PS S ) ) 

L R Tv F-rFFF / 4 • 

IMfcFMi.CE.tHRMN) CCTO 250 



TE IS F C I N T EAS TEE SMALLEST CRRCR FCUNO SC FAR, SC IT 
ANT REFLATES TFF PREVIOUSLY FOUND BEST PCINT 



IS SAVED 



PS , PT = TFt PEST STATIC# TCTAL FRCSSURF FECND 

XM IN , VMK •= TEE Y A to » PITCE ANCLES WHERE THE MIN I N L N C F R C R to. AS FCUNC 

J..>MN = EFn, 

P s= F s S 
‘ - T = r- T T 

XM I E= > 

YM I E=Y 

5 C X-X-iDEL 

IE ( >#LE#ANAX) CCTC 170 



00000950 
30300970 
0 00 00930 
0 0 3 C 0 9 9 0 
30301000 
930C1Q10 
CCOC1 020 

: c o c i o 3 o 

30)01040 
00001 050 
0 C O C 1 C6 0 
OOOCl 070 

: co c i cao 

30001 090 
•3 0 0 C 1 10 0 
0 0 3 01 1 1 0 
OOOCl 120 
0 03 0 1 1 30 
0300 1 140 
3 0 0 0 1 1 5 0 
0 0 0 0 1 16 0 
CCOC 1170 
OOOCl 130 
3 0 0 C 1 19 0 
3 C 0 C 1 200 
COO C l 2 1 0 
00301 220 
003 0 1 225 
00001230 
0 CO 01 24 0 
30001 250 
00001 260 
COCCI 270 
00001 230 
COO 0 1 290 
30001300 
00301310 
00 3 0 I 32 0 
CCOC 1 330 
0000 l 340 
OOOCl 350 
COO Cl 360 
COO C 1 370 
00001390 
C CO C 1 290 
30301400 
00001410 
00001420 
0000 1 430 
COOOl 440 
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err. n r. n r r > n r ( ' r c r r r r r, r r r r r r 



VELOCITY 



300 Y=Y*CEL 

IF(T.Lf:.F*A>) C( TC 1 6 C 

Ati s.CM I Mf REDUCING THE GF 10 SIZE uMIl ThC f F f* t F IN THE 

ANCLl flaots c.l degrees 

IMCCL.LE.C.5CI) C C T C 3 f 0 

'At '^lPlAI THE RRCCCCuRL ARCLMJ TIE ElS T PCINT FCLNC SC FAR 
EXClPT USING a CR1D 1/3 AS h*DL 

AN I h = XN IN-CCL 
AM>->fMN i DLL 
F M l I = > VI N - OF L 
PNA>=YNIN + CEL 
DEL = OCL/2. 

GLTC 1 5 C 



CAUU ATC IFF CFEIRFD CJAMITIES. FIRST O-FCK INC F C R TFESC CRIERS 
lt-L A u NFANS TFE RANGE CF THE CALIJRATICN TALLE WAS EXCEEDL9 
II TFE LAST SCAN 

STATIC PFFLSURl <= 0, TK KG 1C VELOCITY RECUIELS A PCSITIVF 

STATIC PRESSURE 



PF C = 
VE L = 
GC 

>A A C F = 



f L L I C CCNGITV (KG/MM2) 
fLUC VELOCITY {V/ SCC) 

‘cmc \nccnr cr fli io <n/scc) 

FLUID M(h IUVElR 



330 UMFL.NE.C) VFITE<G,7CC0) 

/ZCC i-LKAT(<-XM WARNINC THE RANCL CF THE CALIERATICN 
1 MALLE MCHT NET FAVC CEEN SJD IC IEM TC M 

; MLLCW F f C PE f CALCLL A T I CNS M 

I H c £ «LE « C • ) GCTC 450 
I, EC = P£ / <RC4£KLVPMK4^72.U;> ) 

V AC h = SQP T ( < (PT / PS ) **< (GANNA— I • ) /CA^MA )- l •)/( (GANNA-I . )/2 . ) ) 
CO = SCRT < C ANNA* RCAS*(TC*273. id) 
v LL =C C * N A C I 

aRITEIfc.EClC) F S • PT • X N I N # Y N I N • V FL t N A C F 
6C I 0 PCh I A T < IX ,ZF 1 2*2.ex,2Fe.2*SX,Fe.2.F3.3> 

GCTC 1C 



A ILLATIVE STATIC PRESSURE FAS dCEN FCUNC 



430 *K*1£(£'7CIC) r£, PT.XNIN.YNIN 

7010 F C R N A T ( • ICCATIVE STATIC P N CS SUR E • • / . 

1 * FS.Pl .YA^PITCF IM4F12.2) 

C.U 10 

ICCC t. R I I E ( N , 7 C 3 C ) 

7030 f-CKNAU • AN INPUT ERFvCR CCCLRREO WHILE • • 

1 •RfACiNC IN TFL PRCCE CHARACTERISTICS’ ) 

999 S FC F 
t NC 



COO 0 I 450 
C CC C l 46C 
C C 0 0 i 4 7 3 
CCCCl 480 
w. 0 301 4 9 0 
C C J 0 I S 0 0 
C C 3 C I 5 I 0 
C 3 ) 3 I 320 
C C 3 C I 5 3 0 
C C 3 0 l 5 4 0 
C C 3 C I s 5 3 
C 0 3 0 1 5 5 0 
CC0CI570 
C 3 0 C l 5 8 0 
C C 3 C I 3 9 3 
C C : C l CO 0 
CC 3 Cl 6 10 
CCC C l C2C 
CCOCICJO 
C 3 ) 3 1 E 4 0 
CCCCl £45 
C C 0 C I e 5 J 
COOOiC bO 
CC301C7C 
CCCOIC^C 
C 00 01 730 
C 0 0 3 I 7 13 

c c : e i 7 2 c 

CCO 0 I 730 
C CO C I 74 0 
C C 3 C l 7 5 0 
C0J01 7o 0 
CCCCl 7 73 
C 00 C I 780 
C 03 0 I 790 
COCCI see 

003 ci e i o 

C030I 820 
CC0CI83C 
COCCI 840 
CCO Cl 6 5 C 
C 0 0 0 I 8 6 0 
C 30 0 l 90 0 
C CC Cl 91 C 
L 00 C I 920 
C 30 0 I 9 30 
0 30 01 940 
C03C195Q 
00 301 960 
CCO 01 9 70 
3C3C1960 
00301990 
CC002C 10 
C 00 02 02 0 
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r>< r non r> n r n r rnmnc nr, 



velocity 



’2UCk‘ L 'TINC I NPT 

TP ’ 1 i Hi3i: r .UINf HEADS IN THE DATA fTT< THE rRCEE CP A J A C T Z I '> T I C 2 • 
IT cAN bl. CHAN CL 13 TC ANOTHER jL I T ALL E FCRV IP EIICU IRCE. 



' \ \ # N ( J = 
rlL(M) •’ 

l 3 iv.i ( * ) - 

CP< NA ,NP ) = 



N U 1* E L P CT POINTS GN TPE AL^I-A, Phi AXIS 
ALPHA VALUES ON THL A< IS nr T P i£ L fi L[FRATIC^ 

t a o l f: 

PHI VALUES ON TEE AXIS CP THE CALIOPATICN 
T A £L C 

MATRIX CCNTAINING T h F VALUES CP CP FCP THE 
PARTICULAR PROOF 



d J J 0 

8 0 10 
8C CO 



3U~ ROUTINE 1RPT(NA#NP*PK3.CP. IEP) 



s i v e n c i c n r r h ( i q • 2 ) , 

K L A [ ( l » SC C C . F N C =0 0 9 ) 
PL i' NAT (CIA) 

ut: Acu.ecic. Nc-s^n 

P L A f ( i .PC 10 • C N C = 3 S S ) 
f L R N A T ( 1 c F a . 3 ) 

n e a c ( p . sc 20 •LNc=qqq ) 
r o n v a t ( i o h r . o > 

I E E = c 
L T URN 



C p C 1 9 . 1 9 ) 

N A * NP 

(PRIM I . I J.I-l.NA) 

( p R n ( J . 2 ) . J = 1 * N P ) 

( ( CP( 1 . J ) # J= 1 .UP ) . I 



NA ) 



T I 

SC 



RE PA£ LLPH AN CRRCR" WHILE 

an hirer flat, icr, is set 



INPUTTINC TPE CMA, 
= 1 



RETURN 

l.mo 



C C 0 C2 C 3 C 
00002043 
00002050 
0 C 0 0 2 C 6 0 
0 C 0 0 2 C 7 0 
J0002080 
3COC2C90 
0 0 0 0 2 1 ) 0 
r3C.MI1 
COO C2 1 20 
0 0 3 0 2 I 3 0 
30002 140 
000 02 150 
30002160 
0 C 0 0 2 1 7 0 
00002180 
: 3 o c ? i q o 

: 3 0 C 2 2 3 C 
0 0 3 C 2 2 1 0 
CCCC2220 
30302230 
00002240 
: 00 C 22 UC 
) 0 3 02260 
30002270 
33002280 
30302233 
00 3 02 JO 0 
3 C 3 C 2 3 1 C 
33002320 
303C2JT0 
30302340 
00002350 
30002363 
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